A Flattened Protostellar Envelope in Absorption around L1157 
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Deep Spitzer IRAC images of L1157 reveal many of the details of the outflow and 
the circumstellar environment of this Class protostar. In IRAC band 4, 8 ^m, there is 
a flattened structure seen in absorption against the background emission. The structure 
is perpendicular to the outflow and is extended to a diameter of ~2'. This structure 
is the first clear detection of a flattened circumstellar envelope or pseudo-disk around 
a Class protostar. Such a flattened morphology is an expected outcome for many 
collapse theories that include magnetic fields or rotation. We construct an extinction 
model for a power-law density profile, but we do not constrain the density power-law 
index. 

Subject headings: stars: formation; stars: circumstellar matter; stars: pre-main se- 
quence; infrared: stars 
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1. Introduction 

The L1157 dark cloud in Cep heus (IRAS 20386+6751) conceals a young protostar, a so-called 



( Gueth et al.1 2003 



Beltran et al 



Class source (lAndre et al.lll993 ) , whi ch is deeply embedded within a large circumstellar envelope 



20041 ). LI 157 has a large po werful molecular outflow that is the 



prototype of chemically active outflows (jBachiller et al.ll200ll ). Despite the attention that L1157 
has received at radio wavelengths, few observations have be en made in the near to mid-infrared 



outsi de of observations of the outflow in H2 and K-band (e.g. iDavis &: Eisloeffel 



1995 



Cabrit et al 



19981 ). Only re cently have sensit ive instruments been available to observe these objects shortward 



of 10 /im (e.g.. lTobin et al.ll2007l ). The outflow carves cavities in the circumstellar envelope, which 



allow photons from the embedded central source to escape and scatter off dust in the cavity at NIR 
wavelengths. The morphology of the scattered light can be used to probe many of the fundament al 



properties of the source such as opening angle, envelope mass , etc. (e.g.. IWhitnev et al 



2003 



Tobin et al 



2007 



Robitaille et al.ll2007l : ISeale fc Loonevi 120081 ). 



In this letter, we present new, deep Spitzer Space Telescope observations of L1157. The IRAC 
continuum emission is dominated by molecular line emission in the outflow. Near the source there 
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is a small amount of emission that may be attributed to scattered light and perhaps molecular line 
emission that is highly excited by the outflow jet. Other than the enormous outflow (~0.5 pc) to 
the north and south, the most prominent feature observed is a large, flattened absorption feature 
at 8.0 /im and less defined at 5.8 (im. This absorption feature is a flattened circumstellar envelope 
observed in silhouette against the Galactic infrared background. 

The distance to L1157 is important to any physical interpretation, but the distance is highly 
uncertain. The molecular clouds in Cepheus have three characteristic distances, 200, 300, 450 pc 
(|Kunlll998l ); L1157 has a similar galactic latitude as the 200 pc and 300 pc absorbing clouds. Due 
to this, we disagree wi th the c urren t accepted distance of 440 pc. This value was based upon a 
study of NGC 7023 in IViottil (|1969i ); L1157 is not in clear association with this cluster. In this 
letter, we use a distance of 250 pc. 



Observations 



The biggest breakthrough in observing Clas s sources in sca ttered light has come with the 
sensitivity of the In frared Array Camera (IRAC) (jFazio et al.ll2004l ) on the Spitzer Space Telescope 
(jWerner et al.ll2004l ). This enables observers to see through the dust enshrouding a Class source 
and reveal structures close to the embedded source via scattered light and emission lines. IRAC 
has channels numbered 1 through 4, corresponding to central wavelengths of 3.6 //m, 4.5 //m, 5.8 
fj,m, and 8.0 ^m, respectively. 

The L1157 cloud was observed with IRAC on 2006 August 13. Observations were carried out 
using the High Dynamic Range mode with frame times of 30 seconds using a cycled dither pattern 
of 30 positions and a small scale factor achieving a total integration time of 900 seconds. This 
observing sche me increases the o verall sensitivity to scattered light in the outflow cavity of Class 
sources (e.g.. iTobin et al.l 120071 ). Post-BCD pipeline products (version S14.4.0) were solely used 



in this letter. A color-composite image, with exaggerated band 4 stretch of the field is shown in 
Figure [TJ In addition, a greyscale image using only IRAC channel 4, 8 ^m band, is shown in Figure 

El 



3. Results 



3.1. Flattened Envelope in Absorption 



Even a cursory examination of Figures U and El rev eals a clear absorpt ion feature perpendicular 



to the outflow. Much like the proplyds in Orion (e.g.. lO'dell et al.l ll993). the structure is seen in 
absorption against a bright background. In our case, the bright and pervasive background emission 
is from the 8.0 fim band. There is also a lower level of intrinsic background emission in the 5.8 /im 
channel, which shows a less prominent absorption feature. However, the level of absorption is too 
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low for a rigorous analysis. If the absorption is due to the opacity of the circumstellar material, 
we would also expect to see absorption features in the 3.6 and 4.5 but in those channels, 
the background radiation field is nearly zero. The most likely explanation is that the background 
emission is from Polyc yclic Aromatic Hy drocarbon (PAH) features, which are strongest in the 5.8 
and 8 /tin bands (e.g., iFazio et al.l 12004 ) . PAHs are typically fluorescing due to the absorption of 
ultraviolet photons from the surrounding interstellar r adiation field; thus, PAH emission is a ubiq- 
uitous feature of the diffuse interstellar medium (e.g., iFlagey et al.l 12006). So while the proplyds 
in Orion are seen against the bright nebula emission, L1157 is seen against PAH emission. The 
L1157 absorption feature is large, ~l-2' or ~15, 000-30, 000 AU at a distance of 250 pc (see JQ) and 
flattened with an axis ratio of ~0.25. Due to the size and lack of any kinematic information, we 
will call it a flattened envelope or a pseudo-disk, not a circumstellar disk; such flattened density 
enhancements are expected from many theoretical constructs that include magnetic fields or ro- 



1993 



Hartmann et al 



1994 



tation in the collapse process (e.g.. iTerebev et al.l Il984l ; Fiedler &: Mouschovias 



1996). 



1993 



Galli fc Shu 



The observed extinction is probably dominated by a combination of silicates and ice. The 
8.0 /im channel overlaps with the 9.7 /mi silicate feature in half of the bandpass. However, ice 
features are present in both the 5.8 and 8.0 /im channels: the 6.85 /mi (CH3OH or NHJ) and the 
9.0 /im NH3 ice features in the 8.0 /im channel, and the 6.0 /im H2O ice feature in the 5.8 /mi 
channel. The qualitative appearance of the absorption does correl ate very well with observations 
of the dust continuu m emission structure s een at A = 1.3 mm (jGueth et al.l 120031 ) . as well as 
the ammonia emission Bachiller et al. (j 19931 ). but as these trace the dense core, this may not be 
especially surprising. 



3.2. Outflow Features 



L1157 has one of the most studied and well-developed bipolar outflows. The outflow has an 
inclina tion of ~80° and a slower blueshifted (southern) lobe than the redshifted (northern) lobe 



(e.g. iGueth et al.l 119961 ; iBachiller et al.l l200ll ) . The S shape morphology and the three peaks at 
point reflection symmetry seen in CO and SiO emission are we l l explained by outflo w precession 
with a cone angle of ~15° (|Gueth et al.lll996l ; Izhang et all I2OO0I ; bachiller et alJboOll ). 



The IRAC data have a remarkable coincidence of emission structure with the molecular out- 
flows, Figure [21 We measure the largest separation of peaks as ~15° and an angle of the outflow 
extension in width ("east-west" direction) a s ^35° . This means that the precession cone has a ~15° 
angle, consistent with previous studies (e.g. lZhang et alj|2000i ; IBachiller et al.ll200ll ). and that each 
episodic shock has ~10° opening angle. The precession cone angle here indicates the total angle of 
the cone, and the opening angle represents half of the outflow opening. From the peak positions, the 
precession period is estimated as ^3050 years, ass u ming 250 pc distance and 100 km s _1 constant 



outflow velocity from the model of IBachiller et al.l (|200ll ). 
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Unlike the extinction structure, the outflow features are commonly shown in all four IRAC 
bands. The excitation mechanism for the broad energy range of emission is beyond the scope of this 
letter; it requires spectroscopic observations and detailed modeling to understand level populations 
in a large energy region. However, we can cons ider hydrogen recom bination lines (the bipolar 



outflow regions were dissociatively shocked, e.g., iBachiller et al.ll200ll ). and molecular hydrogen 
ro- vibrational and rotational lines (based on the chemical complexity of the outflow). In fact, 
many of these transitions are in the IRAC bands with some metal ion fine structure lines (e.g. 



Rosenthal et al.1 12000 : iNoriega-Crespo et al. 



2004; 



Neufeld et al.ll2006l ). The 3.6 /xm band may be 



dominated by H 2 1-0 0(5,6,7) and 0-0 S(13). IRAC bands 2, 3, and 4 may be dominated by H 2 
0-0 S(12,ll,10,9), S(8,7,6), and S(5,4), respectively. Some hydrogen recombination lines such as 
Brackett a in band 2 and ion fine structures such as [Ni II] , [Ar II] , and [Ar III] in band 4 also may 
contribute. Finally, PAH emission contributes in bands 3 and 4. 



3.3. Optical Depth Measurement and Models 

The absorption feature of the flattened envelope structure in L1157 provides an excellent 
opportunity to examine the envelope material by measuring the optical depth along the feature. 
This is especially true as the L1157 flattened envelope is nearly edge-on, isolated, and surrounded 
by relatively smooth background emission. The background and associated average uncertainty was 
measured using the IRAEj task "imstat" in three areas within 90" from the central source (0.54 ± 
0.04 MJy/sr). The areas were free of stars and not significantly affected by residual absorption from 
the outer envelope. A constant background may not be entirely realistic, but we di d explore the 



possib ility of modeling the background using median filtering in a similar method to I Simon et al 



(|2006l ). However, using a background model did not improve the analysis. 



The intensity of the absorption feature was measured perpendicular to the outflow, radially 
away from the central source. This was done using SAOImage DS9 to measure a "projection" of 
82 pixels in length and 3 pixels in width. The intensity was taken as the average of the 3 pixels (in 
width) for a particular position along the projection. To determine the intensity uncertainty, we 
used the pixel values provided by the uncertainty frame from the Spitze r pipeline. Note, that the 



point spread function of IRAC at 8.0/im is l'/9, and the pixel size is l'/2 (jFazio et al.l 120041 ) 



In order to quantitatively compare the absorption to axisymmetric models of circumstellar 
envelopes, we created a simple model of the opacity through an absorption slice, perpendicular to 
the outflow. The model consists of an edge-on flattened object, i.e. disk-like, with a radial density 
profile, p(r) = /0o( r / r o)~ p 5 where ro and po are the radius and the density one pixel away from 
the center. Using a Cartesian convention where location along the absorption feature is x and the 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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distance along the line of sight is y, one can re-write the density profile along the absorption as 
p(x,y) = po{(x 2 + y 2 )/Vo)~ P//2 The opacity is then of the form 



WR 2 -x 2 

t(x) = 2*wo / ((x 2 +y 2 )/r 2 )-P/ 2 d(y/r ), (1) 



o 

where R is the outer radius of the disk and k\ is the dust opacity. This is compared to the observed 
opacity of r(x) = — Zn(intensity(x)/background) in the average along the projection. This is an 
edge-on approximation; the source has an ~80° inclination, but as we are averaging over three 
vertical pixels, the effect is minimized. 

To fit this model to the data, a grid of models was used for the following variables: the 
background level (although we measure a background level, we still use a grid of values around the 
measured value), a constant for k\ and po combined, the outer disk radius (R), the location of the 
structure center in pixels (Xq), and the power law dependence (p). Note that there are no a priori 
assumptions as to the values of n\ or p$. This parameter space is then compared to the 8.0 pm 
data using a x 2 likelihood. The pixels close to the central source were not used in the fits, as they 
contain bright emission from the scattered light close to the protostar, see Figure El 

Our model parameter grid used background levels = 0.45 to 0.60 MJy/sr in steps of 0.01, 
p=0.5 to 3 in steps of 0.5, R= 18" to 84" in steps of 1.2", X was fixed to 0" or ±1.2" of the peak 



of the compact dust emission (IBeltran et al.ll2004l ). and the constant=0.01 to 2.0 in steps of 0.01. 



The model had nearly 3 million grid points. Using these parameters, the simple model of optical 
depth successfully fit the data with high confidence levels. Models were considered good fits if 
the confidence level was >90%, i.e. >10% likelihood (see Table H]). In general, the fits are not 
well constrained. The power-law index, p, fits range from 0.5 to 2.0, and the outer radius, R, was 
constrained to >27.6"; however, the maximum is not well constrained as we only modeled ±50" of 
data. Examples of the "best fits" for each acceptable density power-law are shown on the 8.0 pm 
data in Figure [3l 



4. Discussion 



The shape of the absorption feature is especially intriguing, as it looks like a disk structure 
perpendicular to the out flow axis. Th i s is the first clear detection of a flattened envelope or pseudo- 
disk in a Class object. Galli <fe Shu (jl993al fbl) have shown that a modest magnetic field structure 
modifies infall from the initial spherical cloud to form a so-called "pseudo-disk" ; a flat thin structure 
in the equatorial plane that is not rotationally-supported, thus collapsing. This type of structure 



is also seen from simple flat sheet models of collapse (e.g. jHa rtman net al.l [1994J, 119961 ) , as well as 



detailed ambipolar diffusion models (e.g.. iFiedler fc Mouschovias!ll993l ). On the other hand, this 
structure is large ~15,000-30,000 AU, depending on the background used. That size is somewhat 
larger than the inn er envelope size estimated from interfero metric models of th e dust continuum 
(jLooney et al.ll2003l ). However, the single-dish dust emission (jGueth et al.ll2003l ) is extended along 
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the same axis as the absorption, which argues that the inner envelope in the equatorial plane either 
has higher density, and/or different dust opacity properties. 

Our modeling results show that the properties of the structure, as determined by the absorption 
model, are consistent with the above theoretical constructs, i.e. flattened envelopes and density 
profiles. Although we model a range of indexes ranging from p = 0.5 to 3, only the 0.5 to 2 provide 
acceptable fits at the 90% confidence level with the vast majority of fits being p=1.5. 

To better explore t he physical meaning the model, we assume dust opacities (dust plus gas) 
from lLi fc Draind (hoOll ) of i^s.Ofim — 5.912 cm 2 g 1 Although using interstellar dust opacities for 



, — — , „ 

Loonev et al. 


2003; 


Natta et al. 


2007) 



assumed ns.o^m, the derived range for the density reference, po, or 1 pixel (l'.'2) from the center of 
the envelope (i.e. 300 AU at 250 pc) is listed in Table [TJ 

In addition, we can estimate the absorbing mass of the flattened envelope component for each 
"best fit" model of Figure [3] and a height of 3 pixels (the size of the box we averaged over). We 
assume that the vertical density profile is constant for the mass estimate, even though the observed 
absorption falls off vertically with scale heights of ~3-4 pixels, using a Gaussian vertical structure. 
Our mass estimate, listed in TableHJ ranges from 0.08 to 0.16 M Q . Without using any model, we can 
also estimate the mass required for the observed extinction. We use the above value for Ks.o/im and 
a background of 0.54 MJy/sr to calculate the mass necessary for the extinction of all pixels in the 
central region below 0.458 MJy/sr (577 pixels). The total absorbing mass required for those pixels 
is 0.19 Mq. It is important to note that this mass is only in the absorbing pseudo-disk, but the mass 
is comparabl e to th e 140" extended envelope detected in the millimeter continuum as discussed in 
Gueth et al.l (|2003l ) with an estimated mass of 0.7 M Q at an assumed distance of 250 pc. This 
implies that although a large fraction of the mass of the envelope is in the flattened structure, most 
of the mass is more diffuse. Due to the log nature of the mass absorption, the density contrast from 
the center of the absorption feature to slightly offset from the absorption feature is a pproximately 
an ord er of magnitude. This is more than expected from th e numerical models of 



Galli fc Shu 



()1993bl ) , but somewhat con sistent with the ambipo lar models of lFiedler k, Mouschoviad (|1993l ) and 
the sheet collapse model of lHartmann et al.l (|1996l ). 



As can be seen in Tabled] and Figure these data can not constrain the model density profiles 
directly; multiple power-laws are allowed. On the other hand, although many theoretical models 
do suggest a flattened envelope structure, the power-law index in the pseudo-disk is expected to 
evolve. In that case, we would not expect a single power-law to well describe the envelope density. 
In fact, the ambipolar models (e.g. Tassis Mouschovias 2005al lbl) suggest that the power-law can 
be episodic in the flattened envelope. Further studies with increased sensitivity need to be compared 
directly to the theoretical density profiles in the flattened objects to say anything more completely. 



We thank F. Gueth for providing the A = 1.3 mm continuum emission map and R. Bachiller for 
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providing the CO 2-1 emission map. This work is based on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology 
under a contract with NASA. Support for this work was provided by NASA. This research has 
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Table 1. Summary of Acceptable Fits 



Power law 


Background 


R 


Constant 


Po a 


Mass b 


Max. 


Number 


P 


(MJy/sr) 


(") 


(kaPOH)) 


(10 -18 g/cm 3 ) 


(M ) 


prob. (%) 


of fits 


0.5 


0.47 - 0.60 


27.6 - 50.4 


0.02 - 0.03 


0.75 - 1.13 


0.08 


99.2 


130 


1.0 


0.46 - 0.60 


28.8 - 84.0 


0.06 - 0.13 


2.26 - 4.90 


0.16 


99.9 


1820 


1.5 


0.46 - 0.60 


30.0 - 84.0 


0.26 - 0.56 


9.80 - 21.11 


0.14 


99.9 


7253 


2.0 


0.47 - 0.51 


34.8 - 84.0 


1.04 - 1.63 


39.20 - 61.43 


0.15 


48.3 


2622 



a Assuming d=250 pc, for other distances multiply by (250/d). 

b The mass is estimated using the "best-fit" models from FigurcQ5]with d=250 pc. For other distances multiply 
by (d/250) 3 . 
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Fig. 1. — Color IRAC image of the L1157 region with Chl-blue, Ch2-green, and Ch4-red. The color 
stretch is slightly exaggerated to emphasize Ch4 (8 ^m band) where the extinction is the largest. 
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Fig- 2.- 



IRA C ba nd 4 greyscale image of L1157 overlaid with the CO 2-1 emission from 
Bachiller et al.l (120011 ). The upper-left inset is the same greyscale without the contours. The 
lower-right inset is a closeup (±75") of the absorption structure with contours of 0.042 MJy/sr x 
1, 2, 3, and 4. 
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Fig. 3. — "Best fit" examples for the 4 acceptable density power-law indexes on the 8 fim data. 
The fit parameters are: (p=0.5) background = 0.5 MJy/sr, constant = 0.02, Xq = -1.2", and R 
= 39.6"; (p=l) background = 0.57 MJy/sr, constant = 0.11, X = -1.2", and R = 55.2"; (p=1.5) 
background = 0.49 MJy/sr, constant = 0.37, X = -1.2", and R = 55.2"; (p=2.0) background = 
0.48 MJy/sr, constant = 1.3, X = -1.2", and R = 40.8". 



